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Abstract 

We give a complete numerical description of the geometry of the four-point contact inter- 
action of closed bosonic string field theory. Namely, we compute the boundary of the relevant 
region of the moduli space of the four-punctured spheres, and everywhere in this region we 
give the local coordinates around each punctures in terms of a Strebel quadratic differential 
and mapping radii. The numerical methods are explained in details. And the results are 
translated into fits, which can in principle be used to compute the contact interaction of any 
four off-shell string states. 
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1 Introduction 

Although a field theory of closed bosonic strings is known (PUEJIH]), until now only very little 
progress has been made in understanding the fate of the bulk tachyon of the closed bosonic 
string theory. Yet, calculating the tachyon potential, and rewriting the theory around one of 
its minimum, is probably the first thing that one should do if one is given a field theory with 
a tachyonic instability. The problem with closed bosonic string field theory (CSFT) is its 
complexity; unlike open string field theory, which is cubic in the string field, CSFT is non- 
polynomial. This adds a great deal of complexity, because all terms higher than cubic will 
involve an integration over a subset of the moduli space of the relevant Riemann surfaces; 
and at every point of this restricted moduli space we need to specify local coordinates, so 
that we can properly insert the off-shell vertex operators. 

It may be nevertheless feasible to do calculations in CSFT, by truncating the action to 
some finite polynomial order; it is possible that the results computed from the truncated 
action converge as the truncation order is increased. If we can do so, one of the most 
important questions that we would like to answer, is whether the tachyon potential has a 
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minimum or not. If it has a global minimum, then it is natural to expect that the theory 
expanded around this minimum may be a superstring theory. The mechanism by which the 
bosonic theory would flow to a superstring theory by tachyon condensation, would have to 
involve spontaneous breaking of the twenty-six dimensional Poincare symmetry. Evidence 
that this can happen in CSFT at low order, was given in [3]. This is also supported by the 
fact that all the ten-dimensional superstring theories can be realized from the closed bosonic 
string theory after some suitable compactification and truncation of the spectrum But 
although the tachyon potential has a minimum to third order in the string field (as was 
originally found in [6J), to quartic order it doesn't have a minimum anymore, as was shown 
in 0. Clearly, to get some new insight into this question, we need to express the tachyon 
potential to higher orders. 

Another interesting question, which has appeared more recently in the literature (see Bj, 
jH] and references therein), is that of the condensation of closed tachyons localized on the 
apex of an orbifold. It has been conjectured that the condensation of such tachyons changes 
the background geometry to that of an orbifold with a lesser deficit angle, until it eventually 
condenses to flat space. A numerical check of this conjecture was attempted in j^j in closed 
bosonic string theory. At low order and truncation level, a relatively good agreement (of 
roughly 70%) was found. Here too, a more thorough study of the problem is very likely to 
require contact interactions of higher orders. 

In [Zj, Belopolsky calculated the four-tachyon contact term. His method relied on analyt- 
ical properties of the quadratic differentials on a four-punctured sphere (although numerical 
computations were still needed in order to extract the amplitude). Unfortunately, these 
analytical methods don't seem to generalize obviously to the contact interactions of more 
than four strings. The goal of the present paper is to provide an efficient numerical method 
for calculating the geometry of the contact interactions, that does generalize to higher order 
interactions. We also give our numerical results, that describe completely the geometry of 
the interaction everywhere in the moduli space. They can thus be used to calculate the 
contact amplitude of any four off-shell string states. 

The structure of the paper is as follows: First, in Section El we briefly review how to 
calculate the contact amplitudes with the help of quadratic differentials, and how the local 
coordinates at the punctures can be computed from the Strebel quadratic differential. We 
review the moduli space of the contact interaction in Section El In section 01 we expose in 
details our numerical method to find the Strebel differential at a given point in the moduli 
space; then we go on and describe, in Section our numerical method to compute the 
mapping radii of the Strebel differential. We then present our numerical results in Section 
13 in the form of fits. And in Section [7J we show that the result of (7j for the four-tachyon 
contact term is confirmed by our results. Finally, Section El is devoted to discussions and 
conclusions. 
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2 The geometry of the contact interaction 



In this section, we summarize very briefly the notions that are necessary to compute off-shell 
contact amplitudes. More detailed explanations can be found in EH HI] . 

To calculate an off-shell contact amplitude of iV string states on an iV-punctured sphere, 
one must insert the vertex operators at the punctures, transforming them according to 
the local coordinates at the punctures. One then calculates their correlation function after 
inserting antighosts. At last, one has to integrate the correlators on the relevant subset 1 Vo,at 
of the moduli space of A^-punctured spheres. The geometry of the punctured sphere that 
determines the local coordinates at the punctures, is given by the metric of minimal area 2 
with the constraint that every nontrivial closed curve has length > 2ir. It follows that the 
geometry of the punctured sphere is that of A" semi-infinite cylinders of circumference 2tt, 
joining along the edges of a polyhedron whose vertices always join three edges. The constraint 
on the nontrivial closed curves ensures that we do not over-count Feynman diagrams. The 
solution to this minimal area metric is given by a Strebel quadratic differential, which we 
now briefly introduce. 

A quadratic differential <p is a geometrical object that lives on a Riemann surface. In 
the local coordinate z, it takes the form tp = (p(z)(dz) 2 ; and it is invariant under conformal 
changes of local coordinates, in the sense that, if w is another coordinate, we have 

<j){z){dz) 2 = <f){w){dw) 2 . (2.1) 

A quadratic differential defines a metric, whose length element is ds = \/\<f>(z)\ \dz\. In 
particular, in a neighborhood of a second order pole z = p, 

<!>(*) = j^ji + 0{{z-p)- 1 ) , (2.2) 

and if r < 0, the geometry is that of a semi-infinite cylinder of circumference — 2nr. The 
quadratic differentials that we will consider here are regular everywhere except at the punc- 
tures, where they have second order poles. We will call the quantity r in (|2.2|) . the residue 
of ip at the puncture p. It is easy to see from f)2.1j) . that the notions of pole and of residue, 
as well as the zeros of tp, are conformal invariants. If <p has a second order pole with residue 
r at z = p, then there is a local coordinate w, defined in a neighborhood of the puncture, in 
which 

V =- 2 (dw) 2 . (2.3) 

We will call such a local coordinate w, a natural coordinate at the puncture z = p. 

1 In general V Si at is the restricted moduli space of the iV-punctured Riemann surfaces of genus g. Although 
we will consider here only the case g = 0, we will keep this notation. 

2 Actually, since in our case the area is infinite, we have to minimize the reduced area. See |H] for more 
details. 
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A horizontal trajectory of ip is a maximal curve on which <j)(z)(dz) 2 > 0, where dz is 
tangential to the curve. For example, the horizontal trajectories of (|2.H|) are circles centered 
at the puncture w = 0. A critical horizontal trajectory is a horizontal trajectory one of 
whose ends is either a zero or a pole of order one. We note here that, as one can easily 
verify, a zero of order n is the source of n + 2 horizontal trajectories. The critical graph of 
<p is the set formed by all critical horizontal trajectories. 

From now on, we will only consider the case of the 4-punctured spheres, which we will 
describe in the uniformizer coordinate z. The last definition 3 that we need is that of the 
Strebel differential. It is a quadratic differential whose four poles are second order poles 
with residue —1, and whose critical graph has the topology of a tetrahedron. In particular, 
the Strebel differential has four simple zeros, which are the vertices of the tetrahedron, at 
which three critical trajectories (edges) are joining. See figure Q for the critical graph of 
a particular Strebel differential. On this figure we have also labelled the lengths (in the 




Figure 1: The critical graph of the Strebel differential, whose poles are at 0, 1, £, and oo; here we have 
chosen the particular value £ = —0.2 + 1.5 i. It has the topology of a tetrahedron, with a puncture on each 
side. Its zeros are labelled Zi, i = 1,..,4. We have also labelled the edge lengths ti (note that two opposite 
edges have the same length). 

99-metric) of the edges of the tetrahedron. From the residue conditions, it is easy to see that 
two opposite edges (which don't touch each other) have the same length. Also we have that 
the circumference of each face is 2n, in particular t\ + £ 2 + £3 = 2ir, and thus we have only 
two independent lengths, for example £\ and £ 2 . 

3 This definition can be made much more general, but for simplicity, we restrict ourselves here to the 
four-punctured spheres. 
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Using the PSL(2,C) automorphism group of the Riemann sphere, we can always place 
three punctures at, respectively, z — 0, 1, and oo. The fourth puncture will be called £. The 
residue condition for the Strebel differential, is satisfied by the general form 

^ _ -(^ 2 -0 2 a 



zHz-iy(z-o 2 z(z-i)(z-o 

-z 4 + a z 3 + (2f - (1 + Oo) 2 2 + a£z - £ 2 
z 2 {z- l) 2 {z -O 2 



(2.4) 



where a is any complex number. The puncture conditions at 0, 1 and £ are straightforwardly 
verified from the above form. For the puncture at infinity, we need to change the coordinate 
to w = 1/z, and we should verify that 4>(w) = —l/w 2 + 0{w~ v ). This is easily done: 



</>(z) [dzf = <p{w)(dw^ 

=► = fa = l/w) 2 = <j>(z= l/w) — 4 = -— 2 + Oiur 1 ) , 

\dw J w A 

where by 4>(z = l/w) we mean the function <f>(z) in the coordinate z, evaluated at z = l/w. 
And we see that (J2.4j) is therefore the most general form satisfying the residue condition. 

The metric defined from the Strebel differential is a metric of minimal area. Moreover, 
by a theorem of Strebel ([H]), the Strebel differential for a given £ is unique. Therefore, if 
all nontrivial closed curves on the punctured sphere, have length > 2tt, we will have solved 
our minimal area problem. One of the results of this paper is to determine numerically the 
coefficient a in (|2.4jl . and therefore (p, for all values of £ in the restricted moduli space. 



The Strebel differential defines (up to unimportant phases) the conformal maps from the 
four punctured disks (with local natural coordinates w n ) to a 4-punctured sphere, in such 
a way that the punctures w n = are mapped to the punctures z = p n , and the boundaries 
\w n \ = 1 of the unit disks are mapped to the part of the critical graph surrounding p n (see 
Figure EJ)- Our notation is that the punctures are labelled by the index n running from 1 
to 4, with pi — 0, p 2 — 1, Ps — £ and p 4 = oo. The method for expressing explicitly these 
conformal maps from the knowledge of tp was described in Tol . First, we expand ip in a 
Laurent series around the puncture p n . 

P= (—, 1 —; + ^^ + b + b 1 (z-p n ) + ..) (dz) 2 . (2.5) 

V [Z-PnY Z-p n J 

The coefficients \ are calculated straightforwardly from ([2.4)1 . The map from the punctured 
unit disk to the punctured sphere, is expanded as follows 4 

Z = K(w n ) =p n + (p„W n ) + d\ (p„W n ) 2 + d 2 {PnWnf + ... . (2.6) 

4 0ur notation here differs from |1U|. It is convenient to expand in powers of (p n w n ), instead of w n . 
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Here p n = is the mapping radius of the conformal map h n . In the natural coordi- 



nates w n , the quadratic differential has the form 



1 



^{dw n f . 



(2.7) 



w 



n 



Now we use ()2.5|h f)2.6j) and ()2.7|) to calculate recursively the coefficients di from the frj's. 
For the first few coefficients, we get 



Note that the mapping radii cannot be calculated from the above recursive method. But 
they can be calculated, for example, by the formula ()5.1|) derived in [TOj. As a consistency 
check of our results, we have calculated the conformal maps as a series (|2.fij) to order 18, for 
the particular value f = -0.2 + 1.5 i. We find a = 1.12245 + 1.23948 i, and 



In figure |2l we have represented the four conformal maps by mapping circles of radii k/10 
for k = 1, .., 10. 

3 The moduli space 

As already mentioned, a point £ belongs to the moduli space Vo,4 of the contact interaction 
(also called restricted moduli space), only if all nontrivial closed curves on its corresponding 
tetrahedron have length > 2tt. One can easily check that this is equivalent to the conditions 



where the £i are the edge lengths (see Figure HJ). We show Vo,4 in Figure El It is delimited 
by three curves £>i, B 2 and £> 3 on which, respectively £i, £ 2 and £ 3 are equal to n. 

We don't need to calculate the Strebel differential everywhere on this space because 
there are some symmetries. First, there is the obvious symmetry Then there are six 




Pi 

P2 
P3 
Pi 



0.69819 
0.88628 
1.34118 
0.46138. 



ii, k, h < vr 



(3.1) 
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Figure 2: The four conformal maps as series (|2.6|l to order 18, are represented for the particular value 
£ = —0.2 + 1.5 i by mapping circles of radii fc/10 for k = 1, .., 10. Also plotted is the critical graph of the 
quadratic differential. 



PSL(2, C) transformations that permute the three fixed poles {0, 1, oo}. They are generated 



In total, Vo,4 is thus made of twelve parts separated by the real axis, the line Re(£) = 1/2, 
and two circles of unit radius centered on and on 1 respectively. In our computations, 
we will therefore only consider the subspace A, which is bounded from the right by {£ : 
Re(£) = 1/2}, from below by the unit circle centered on 0, and from the left and the top by 
the curve B\ (see Figure EJ). 

To find the parameter a in the other regions of Vo,4, we apply, for PSL(2, C) transforma- 
tions £ — > g(£), the transformation law ()2.1|) to <p given in ()2.4j) : and when £ — ► £*, we have 
a — » a*. The transformation law for the mapping radius p n at the puncture p n , is 



by 



and 



91 ■ 

92 : 



(3.2) 
(3.3) 



Pn 



g'(Pn)\ Pn- 



(3.4) 
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Figure 3: The moduli space Vo,4- It is bounded by the three curves Bi, B 2 and B 3 . We need only to 
concentrate on the space A (shaded in darker grey), bounded by B±, the unit circle at the origin, and the 
line Rez = 1/2; all the other regions can be reached from A by PSL(2, C) transformations and complex 
conjugations. 

And it is invariant under complex conjugation £ —>■ £*. Note that if g in (|3.4|) . sends the finite 
puncture p n to infinity, one must calculate its transformation in the coordinate w = 1/z. In 
particular, if g(£) = l/£, the mapping radii at the punctures z = and z = oo are invariant. 

Note that we don't know the analytic expressions of the curves Bi. We thus have to 
calculate one of them, say Bi, numerically. This is done by finding the points where l\ = it; 
the result will be presented in Section |H1 But we can already understand a few points of 
this moduli space. First, when £ = 1/2, we can show that the Strebel differential is given 
by a = 2, and then one can verify that £2 = £3 = it. This shows that at the point £ = 1/2, 
the two curves B2 and B% intersect. And therefore, from the symmetry of Vo,4, one also has 
that, at £ = — 1, Bi and B2 intersect, while at £ = 2, it is B\ and B3 that intersect. The 
knowledge that the point £ = — 1 belongs to B\, will be used in Section EJ when we express 
B\ as a fit. Note that all these points on the real axis are singular, in the sense that the 
Strebel differentials at those points have double zeros, which means that the tetrahedron is 
singular, two of its edges having collapsed to zero length. 

1 Fk 

The point Q = | + ^ z is also interesting. Indeed, when C, = Q, the three punctures 0, 1 
and £ are placed at the corners of an equilateral triangle. In fact, this is the most symmetric 



9 



configuration, whose corresponding tetrahedron is a regular tetrahedron. One can check that 
for this configuration, a = 2 + i. 

4 Computation of the Strebel differential 

Here we explain in details how to find the parameter a in (12. 4j) such that the quadratic 
differential is the Strebel differential for the specified position of the fourth puncture £. 
Once a is found, the mapping radii can be computed relatively easily; this will be described 
in the next section. The numerical results for a and the mapping radii will be presented in 
Section |H1 

We define the complex length of a curve C by £(C) = J c ^P(z)dz or, if C is a segment 
from Zi to Zj, 



Zi ,z 3 ) = / ^{z)dz. (4.1) 



The square root must be chosen such that its branch cut doesn't intersect the path of 
integration. Now suppose that Zi and Zj are zeros of p, and that p is the Strebel differential. 
This means that there is a critical trajectory, that we call C, joining Zi and Zj. From the 
definition of a horizontal trajectory, we have that a/ 4>(z)dz is real on C, and therefore £{C) is 
real (and is equal to the length of C in the ^-metric) . By deforming the path of integration, 
we then have that £(z{, Zj) is real. Indeed, even if p has a pole in between C and the segment 
[zi,Zj], the residue of this pole would change the integral only by a real quantity, because 
the residues (in the usual sense) of \/(f)(z) are all purely imaginary. 

Therefore, in order to find the Strebel differential, we will look for the quadratic dif- 
ferential (p that satisfies the condition that £(zi, Zj) must be real for any two zeros z% and 
Zj. Since the quadratic differentials (|2.4jl that we are considering, are already constrained 
to have residues —1 at the four punctures, we easily see that the above reality condition is 
satisfied if £(zi, Zj) is real for two independent segments; for example 

lm£(zx,z 2 ) = 

1 ' ; (4.2) 
1m£(z u zs) = 0, 

where Z\, z 2 and z 3 are three zeros of ip. Solving these two real equations will determine the 
complex parameter a. 

To solve the system of equations ([4.2)1 . we use the Newton method in two dimensions. 
Namely, we have two real equations for the two real unknowns Re a and Ima. In general, 
the Newton method for the system of n equations of n unknowns f(x) = consists of starting 
with a Xq not too far from the solution, and repeating the iteration 
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until Xi converges. Here J(xi) = [%) is the Jacobian. In order to calculate it in our 

V / X=X i 

situation, we note that for a holomorphic function f(z) we have dl ^}^ = ^- m f'(z), and 
^IHmz = ^ e f'( z )- Our Jacobian is thus 

\ ^fe^ "-^tf 1 J V Im ^^3) Re )> 

and we will need therefore to evaluate -£;£(zi, Zj) as well 

But before we go on and solve the system of equations ()4.2|) . we see that there is a 
problem with the computation of the complex lengths ()4.1|) . namely we need a square root 
whose branch cut we shall never intersect, although we have to evaluate it along various 
paths whose positions we don't know a priori. The solution to this problem, a square root 
without a branch cut, is a central tool of our numerical methods; it is presented in the next 
subsection. 

4.1 A square root without a branch cut 

First we note that in all our calculations that present a potential branch cut problem, 
the square root is to be evaluated inside an integral along some path. It has therefore to 
be evaluated at points that are "close to each other". It is thus quite natural to define 
a continuous square root (that we denote by ^/~) along the path, in the following way: 
It has to remember the last result ^/z n -\ that it calculated; then to calculate the next 
value, it computes the usual square root ^/z^, and chooses its sign ^fz^ = so that 

| ^fz^ — ^/z n -\\ is minimal. 

We illustrate the result of evaluating the continuous square root along a path in Figure 
HI where the path is a closed curve winding twice around the origin; its image is a closed 
path that winds once around the origin. We see that the gain of avoiding the branch cut 
(and therefore defining a continuous function), is at the expense of the sign of the function 
not being globally defined 5 . Indeed, in our example we have Wz\ 7^ ^fz~2 although z\ = 
Practically, the continuous square root is implemented by the following (schematized) code: 

complex (complex z) 
{ 

static complex p = 0; 
complex w; 

w = ^/z; 

if (Re(w) Re(p) + lm(w) lm(p) < 0) 
w = —w; 

p = w; 
return w; 



5 The domain of w/~ should really be a double cover of C — {0}; it would then be a well-defined holomorphic 
function. 
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Fi gure 4: The image by the continuous square root of a continuous closed path winding twice around the 
origin is a continuous closed path winding once. Note that the point Z\ = z 2 that is met twice by the path, 
is mapped to two different values. 

The variable p is declared "static", which means that its value is retained between two 
calls of the function. The square root w is compared with the previous result p. If their 
"scalar product" Re(w) Re(p) + Im(w) lm(p) is negative, we suspect that the branch cut was 
crossed, and thus we flip the sign of w. 

We emphasize that the continuous square root makes sense only when it is evaluated 
along a path that does not meet the origin (practically, it has to be evaluated at successive 
values that are close to each other compared to their distance to the origin). We will thus 
only use it inside an integral and won't bother about the sign in front of the integral. Luckily, 
for the problems we need to solve, like (|4.2J) . we will always be able to go around the fact 
that the global sign is ambiguous. 

Finally, if the continuous square root is evaluated along a non-self-intersecting path that 
does not go through the origin, we can always draw a branch cut that does not meet the path. 
The holomorphic square root with this given branch cut is then equal (up to a global sign) 
to the continuous square root. The latter is therefore holomorphic in some neighborhood of 
the path. 
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4.2 Numerical evaluation of £(zi, Zj) 

Without loss of generality, we take Zi = z\ and Zj = z 2 . We need to numerically 6 evaluate 



Z2 fZ2 



VW)dz = I z(z-l)(z-0 dZ (44) 



where Z\,...,z± are the four zeros of <f)(z), and the integration is along a straight path. 
Parameterizing 

z(t) = ^ + ^t = K + ^(t + l), (4.5) 
we can rewrite this integral as 



Zi-zA 2 f 1 ^/(z{t)-z 3 ){z{t) 



where we have written the usual square root y/l — t 2 ; this is fine since (1 — t 2 ) is never 
negative on the integration path. But it is important to note here that we have made an 
arbitrary choice of sign in (|4.fij) : indeed, since the global sign of the continuous square root 
is ambiguous, there is no way to tell which is the right sign in the identity ^fxy = ± ^fx^fy. 
But since we don't care about the global sign of £(zi,z 2 ), what sign we are choosing is not 
important, as long as we do the same choice of sign when computing ^i(zi, Zj) (see (j4.12|) 
and comment thereafter). 

We want to evaluate the integral (J4.6)) numerically. But we observe that, in general, 
one or more poles will be close enough to the path of integration to destroy the accuracy 
of the numerical integration method. To remedy this, we subtract the finite poles from the 
integrand and add back the result of their analytic integrations: 




1 ' } I ' '-1 \z(t)(z(t)-l)(z(t)-0 ^ x z{t)-p n ) 

^anr./^VT^*}, (4.7) 



-! Z{t)-p 

where p% — 0,p 2 — 1 and p 3 = £ are the finite poles, and 



V (Pn ~ Z 3 )(p n ~ Z 4 ) Yl (Pn -Pm) 1 • (4.i 



m— 1 



6 Putting aside the branch cut problem, the integral (|4.4(l can actually be expressed in terms of elliptic 
functions. We still compute it numerically for two reasons. First, the analytic expression is quite complicated, 
and since we have to evaluate it numerically anyway in order to solve (|4.2f) , we wouldn't gain much. Second, 
(|4.4|l is integrable analytically only because the order of the polynomial in the square root is not greater 
than four (see for example and therefore an analytic evaluation of (|4.4|l would not generalize to contact 

interactions of order higher than four. 
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The signs s n = ±1 are chosen in such a way that 



max 

te[-i,i] 



- z 3 )0(t) - 24) s n r T , 



z(t)(z(t) - l)(z(t) - ^z(t)-p 



is minimal. For the four-point contact interaction, we have three finite poles and thus eight 
possible combinations of signs to try. Finally, we need the expression 

1 VT^ 2 r 1 VT^^ z Y + z 2 - 2 Pn 

at = / at , q 



! z(t) -p n z 2 -z 1 J_ l t + q ' ' z 2 - zi 

21 ^-ivr^x/Ti'^ !^!!!^ )■ («) 



ft , 

z 2 - z\ V * V " ^ " \ -sgn(q) if lmg = 
Now that we have subtracted the poles, the integral 



1 



z(t)(z(t) - l)(z(t) - ^ x z{t)- Vn 



can be computed very precisely 7 with a Gauss method with weight \/I — t 2 (Jacobi weight 
function) (see for example t 13J. We use the method of rank 30 which, as we observe, has 
an accuracy of about twenty significant digits for £ = \ + ^i. 



4.3 Numerical evaluation of -^£(zi,Zj) 

Again, without loss of generality, we take Zi = Z\ and Zj = z 2 . We have 



d_ 

da 



—z(zi, z 2 ) 



*2 ^ 

— y<p{z) dz 
- , da 

1 f z ' 2 1 1 
= - / , —, 77 ~,dz (4.11) 

where we have used (12. 4j) . and the fact that 4>(z) vanishes at z\ and z 2 . Now we use again 
the parametrization (|4.5j) . and we find 

d „, x 1 Z" -1 1 dt , 
£( z z ) = — / , 4.12 

^ l 2 A ^(^)- 23 )(^)-z 4 )v / r^ 



7 To avoid cancellation errors (subtracting two quantities that are large and close to each other), we should 
explicitly cancel the pole in H4.10|) . A simplified example is v/t f +1 — ~ = ^ t ^ 1+1 , and the pole at t = has 
been explicitly cancelled. A similar simplification of the expression (|4.1U|) is possible but it gives complicated 
expressions. Instead of doing that, our program always chooses the two independent paths that are the 
furthest away from the poles, so that the cancellation errors are not too large. 
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This integration 8 can be done numerically with a Gauss method with weight 1 / \/l — t 2 . We 
use here the method of rank 30. 

It is important to note, as already mentioned after Equ. (j4.6|) . that we have made the 
same choices of signs in both the calculations of £(zi, z 2 ) and 4^£(zi, z 2 ), namely 



Z 2 ~ Zi y(z{t) - Z 3 )(z(t) - Zj) 

2 z(f)(z(f)-l)(z(f)-fl 



(4.13) 



The relative signs of £(zi,Z2) and -^l(z\, z 2 ) will therefore be consistent if the numerical 
integrations are done one after the other, and if the path of the second integral is the reverse 
path of the first one (see Figure EI). This is why we have written the limits of integration 
in ()4.12|) in the opposite order as in ()4.7|) . In this way, the continuous square root will be 
evaluated at successive points that are close to each other, and we will thus have avoided 
any branch cut problem. 




Z 3 



Figure 5: The arrows show the order in which the numerical integrations should be done. In order for the 
signs of the continuous square root to be consistent, the values used by the numerical integration methods 
should run continuously along the closed path formed by the arrows, starting and ending at z\. 



We observe that if the initial value of a is not too far from the Strebel value, the Newton 
method converges to the Strebel solution. 

5 Computation of the mapping radii 
5.1 The finite poles 

The mapping radius at the puncture p n is given by (|_10j): 

In p n = lim I Im / a/ (fi(z)dz + lne ] , (5.1) 

where z is any point lying on the closest critical trajectory, e is a positive real number, and 
the path of integration is chosen such that it avoids the branch cuts. Here we will take z 
to be a zero, say z\, that lies on the closest critical trajectory, and the path can be chosen 
to be a straight line if we replace the square root by the continuous square root. Moreover, 

8 We are not subtracting the poles here. Indeed, the accuracy that we would gain would only make the 
Newton method slightly faster, without effecting the accuracy of the final result. 
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since the residue of \/(f)(z) at z = p n is purely imaginary, (|5.1|) doesn't change if we replace 
the lower integration limit by p n + ea, where \a\ = 1. Thus 

lnp n = lim(lm / ^/cj)(z)dz + In e ) , a = -r— — —j > (5-2) 

J Pn +ea J \Zl-Vn\ 

and the sign of the continuous square root is chosen in such a way that the limit exists. The 
above expression is clearly very inconvenient for numerical evaluation. But we can remedy 
this by explicitly cancelling the singularity at e = 0. First, we note that 



lirn(lne-( 1 / ^ - dz + 2(ln2 - 1) + In \z x - p n \ ) ) . (5.3) 

~° V. V V Z 1 ~ Pn Jp n +e a Z-Pn 



And thus 



lim 


^Im 










Im 









lnp n = lim(lm/ ( tftfz) - — = ^ - ) dz ) + 2(ln 2 - 1) + In \ Zl - p n \ 
V<Kz) - I V^zA dz + 2(ln2 _ i) + i n \ Zl _ p I . (5.4) 
Thus 

4 . 

Pn = 2l "Pn X 




x exp Im < (z 1 - Pn) 7 7T7 tt Vt dt 

Jo \ Z{z-l){z-i) Z-p n 



where we have parameterized z = Z\ + t(p n — z\). The sign of the continuous square root is 
chosen so that the integrand is finite for z — > p n . The integral is computed numerically with 
a rank 30 Gauss method with weig ht Vt. 

We haven't yet said how to know which zeros lie on the closest critical trajectory. One 
method would be to actually compute the critical trajectory, but this would require solving 
a differential equation, which can be hard numerically, and is certainly time consuming. 
Another method, that we are using here 9 , consists of calculating the radius for all four zeros 
Zi. One of the four results may be different from the other three, we discard it since it must 
correspond to the zero not lying on the closest critical trajectory. Practically, we discard the 
value that is the furthest away from the average value. The redundancy of calculations can 
then be used to evaluate the numerical errors affecting the mapping radii. If we denote by 
Pn the mapping radius for the pole p n calculated with the zero Zi, and if Zj is the discarded 
zero, we identify p n with the average p n ~ p~ n — | X^ =1 Pn ■ And we can estimate the 

variance of this sample a 2 ~ \ X^ =1 (Pn — Pn) 2 - 

12tl 

9 This particular method does not generalize to vertices of higher order. Indeed at order five we already 
have six zeros, and discarding three zeros out of them is thus impossible. For these higher vertices, we can 
use a fixed labelling of the zeros and poles and always make sure that the topology of the critical graph 
doesn't change, so that we always know which zeros are associated with which poles (see |14|). 
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5.2 The pole at infinity 

Here we must calculate the mapping radius in the coordinate w — 1/z, in which the pole is 
at w — 0. From (p(w)(dw) 2 = (f>(z)(dz) 2 we have that <f>(w) = 4>{z = l/w)/w 4 . Thus 




where z\ is a zero that lies on the closest trajectory to the pole at infinity. With the 
parametrization w — l/z\ — t/zi, we find 

4 f r / 1 f l ( ^/-(l-wz 2 )(l-wz 3 )(l-wz 4 ) i 

P± = ~2\ — T ex P Im 1 / T\ vn F\ 

e \zi\ \ [zi J \ w(l-w){l-w£) w 

This expression is numerically calculated in the same way as (J5.5j) . 




6 The results as fits 

The results obtained by our program are lengthy tables of data, namely the parameter a 
and the mapping radii p n given at arbitrarily many points of the portion of moduli space 
A. We do not write here any such table, but instead we express our results in the form of 
fits. Although some precision is lost with the fits, they are reasonably easy to use. The fits 
that we give below all have an error of the order of 0.1%, which should be precise enough 
for many calculations. 



6.1 The boundary of the moduli space 

The element A of moduli space that we are concentrating on, has Re£ < 1/2, Im£ > 0, and 
it is bounded by the unit circle centered at the origin, and by the curve B\ (see Figure EJ) • 
To describe A, we thus need to describe the part of E\ with Re£ < 1/2 and Im£ > 0, which 
we denote by B\. We parameterize this curve in terms of the angle 9 (see Figure with 
the expression 

£(0) = -r(0) e"* , (6.1) 

where 9 G [0, #o]> an d 9 = 1.84644 corresponds to the corner P with Re£ = 1/2. We will 
use the following fit 

£ fit (#) = - r ^(0)e- td , (6.2) 
where r fit (6 l ) is a polynomial of order five in 9 

5 

r flt (0) = l + J2ri9\ (6.3) 

i=i 
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-1 







Figure 6: The parametrization of the curve B\ with the function r(0). 



Note that the leading term has been set to one. In this way, the point £ = — 1, which we 
know is a corner of the moduli space (see Section El) and belongs to Bi, will also belong to 
the fit of B\ (at 6 = 0); and the corner will be fitted properly. 
The coefficients in (|6.3|h are found by minimizing 



j^{9)-r{9)fds 

Here ds is the length element; the denominator in the square root of ()6.4j) is thus the length 
of B\ and has been put there in order to normalize the error. We find 

r &t (9) = 1 + 0.196639 9 + 0.264670 9 2 - 0.130932 9 s + 0.062563 9 4 - 0.015452 9 5 , (6.5) 

with error e r = 0.00022. It is also useful to measure the maximum error made by this fit 

S r = max {\r Rt (9) -r(G)\\ = 0.00088. (6.6) 

6.2 The Strebel differential 

We write £ = x + yi. And we use a polynomial fit of total order 5 in a; and y 

« fit (0= Yl "■r'- i !l J - ( 6 - 7 ) 

i+j<5 
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The coefficients Qj^j £1X6 found by minimizing the error 



J A dxdy 



We find 



a fit 



.3„.2 



(f ) = (0.67867 + 0.46345 i) + (2.07710 + 0.49176 i)x+ (0.20000 + 0.58477i) y 

+ (1.61074 - 0.64246 i) x 2 - (0.09194 - 0.36467 i) xy + (0.05965 + 0.38101 i) y 2 
-(0.33945 + 0.80519 i) x 3 - (0.35448 - 0.94778 i) x 2 y - (1.04747 + 0.85263 i) xy 2 
+ (0.56049 - 0.43233 i) y 3 - (0.63061 + 0.05352 i) x 4 + (0.27490 + 0.76364 i) x 3 y 
-(0.89712 + 0.75355 i) x 2 y 2 + (0.70167 + 0.46374 i) xy 3 - (0.50743 - 0.12128 i) y 
-(0.07990 - 0.08258 i) x b + (0.32821 + 0.08618 i) x 4 y - (0.08918 + 0.20755 i) x 3 y 
+ (0.37843 + 0.20916 i) x 2 y 3 - (0.15046 + 0.08495 i) xy 4 + (0.11992 - 0.00719 z) y 5 

(6.9) 

with e a = 0.0011. The maximum error made by this fit is 

5 a = max{|a fit (0 -a(0|} = 0.0037. (6.10) 

6.3 The mapping radii 

We use again polynomials of order 5 

Pn(0= E Am>*V, ( 6 - n ) 

i+j<5 

but here we minimize the average of the square of the relative errors 



6 Pn 



\ f A dxdy 



(6.12) 



This is more natural than the absolute error, because the mapping radii appear as a funda- 
mentally multiplicative factor fsee 12 .6)) ; in some sense, they measure the sizes of the images 
of the unit disks. The results of the fits are 

= 0.31883 + 0.14609 x + 0.18216 y+ 0.47185 x 2 + O.U573xy + 0.15017?/ 2 
+0.01342 a; 3 - 0.11152 a; 2 ?/ - 0.29992 a;?/ 2 + 0.02631 y 3 - 0.15920 x 4 + 0.01005 x 3 y 
-0.28574xV + 0.17921x?/ 3 - 0.10879?/ 4 - 0.02838 x 5 + 0.08245x 4 y - 0.01874 a:V 
+0.12423 x 2 y 3 - 0.03863 xy 4 + 0.03250 y 5 , (6.13) 
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pj^f) = 0.78098 - 0.59124 x + 0.17253 y - 0.40055 x 2 + 0.44703 xy - 0.10848 y 2 

+0.10574x 3 + 0.27517a; 2 ?/ - 0.03518x?/ 2 - 0.05091?/ 3 + 0.15196 x 4 - 0.08497x 3 ?/ 
+0.05597 xV - 0.08168 xy 3 + 0.07637?/ 4 + 0.01799 x 5 - 0.08381 x 4 y + 0.02112 x 3 ?/ 2 
-0.05276 x 2 y 3 + 0.02554x?/ 4 - 0.01999?/ 5 , (6.14) 

= 0.29797 - 0.36138 x + 0.18030 y + 0.17997a; 2 + 0.49329 xy + 0.32231 y 2 
-0.33706a; 3 + 0.61497a; 2 ?/ - 0.72821a;?/ 2 + 0.23999 y 3 - 0.16258a; 4 + 0.41692 x 3 y 
-0.71139 xV + 0.46090 xy 3 - 0.24944 y 4 + 0.01369 a; 5 + 0.11153 x 4 y - 0.13361 x 3 y 2 
+0.21281 x 2 y 3 - 0.10159 xy 4 + 0.05898 y 5 , (6.15) 

pf(0 = 0.95136 + 0.48365 x - 0.36818 y - 0.16100 a; 2 - 0.52786 xy - 0.01609 y 2 

-0.20922 x 3 + 0.00970 x 2 y + 0.19654 xy 2 + 0.02470?/ 3 - 0.00082 x 4 + 0.17750 x 3 y 
+0.08937x 2 ?/ 2 - 0.0Ul7xy 3 + 0.01641 y 4 + 0.01554 x 5 + 0.00976 x 4 y - 0.03639 x 3 y 2 
-0.03248 xV - 0.00411 xy 4 - 0.00697 y 5 . (6.16) 



The errors are respectively 



e pi = 0.00051 

e P2 = 0.00025 

e P3 = 0.00035 

e P4 = 0.00028 



(6.17) 



while the maximum relative errors 5 Pn = max^^ { Pn*(0 ~~ Pn(0 \ /Pn(0} are 



5 P1 = 0.0020 

& : SJSS <«■»> 

5 P4 = 0.0014. 

7 The four-tachyon contact term 

As en elementary check of our results, we compare the easiest quantity that is calculable 
from our data, namely the four-tachyon contact term, with a result previously obtained, with 
a different method, in It appears as the coefficient V4 in the classical tachyonic potential 

00 N 



N=3 



(note that this is not the effective potential, it is the potential of the tachyon when we ignore 
all other fields). The coefficients are given by ([7j, [TP] ) 



'N-3 \ N 



V N 



JV o,n \n=l / n=l 
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where Vo,at is the restricted moduli space of the N-punctured sphere, and p n is the mapping 
radius at the puncture £ n . Thus 



where we have used the fact that the measure {pip2psp4)~ 2 d 2 £, is invariant under the action 
of PSL(2,C) transformations that permute the fixed poles {0,1, oo}, and under complex 
conjugation; the integration over Vo,4 can thus be replaced by 12 times the integration over 
the twelfth of moduli space A. Using our full data, we get 



in agreement with the result V4 = 72.39 found in jjj. The uncertainty in (|7.4|) was computed 
from the uncertainties in the mapping radii (see Sectional), and the uncertainty from the 
numerical integration over A, which was done with a Monte-Carlo technique. 

It is also interesting to compare this to the same calculation done with the fits given in 
Section IH1 We obtain uf* = 72.41, less than 0.03% away from the value (|7.4|) . 

8 Conclusions 

We have described numerically the geometry of the contact interaction of four off-shell closed 
bosonic string states. Our results give the Strebel differentials and the mapping radii every- 
where in the twelfth of moduli space A, from which we can extend the results to the whole 
moduli space Vo,4 by straightforward PSL(2, C) transformations and complex conjugations. 
We have illustrated how the knowledge of the Strebel differentials and the mapping radii 
can be used to calculate the local coordinates around the punctures of the Riemann sphere. 
Our results can therefore be used to insert off-shell vertex operators on the four-punctured 
sphere, which is needed to compute off-shell contact amplitudes of four closed string states. 

The numerical methods that we have developed can be generalized to contact interactions 
of more than four strings. Actually, the motivation for the present work was to compute the 
bulk tachyon potential and investigate whether it has a minimum or not. While the present 
paper doesn't offer any new information about the tachyon potential, it confirms its quartic 
term, previously calculated in [Z| with a different method. A generalization of the work 
presented here, to quintic order in the string field, is in progress ([14]) and will hopefully 
give the tachyon potential to order five. 

It would be interesting to use the results of this paper for doing concrete calculations in 
CSFT to quartic order in diverse backgrounds, for example in an orbifold background ([HUH])- 
Another interesting application would be to do computations in heterotic string field theory 




(7.3) 



v 4 = 72.390 ± 0.003 , 



(7.4) 



(USD- 
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